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ABSTRACT 


This  paper  deals  with  the  problem  of  selecting  the  best  population  from  among 
fc(>  2)  two-parameter  exponential  populations.  New  selection  procedures  are  pro- 

I 

posed  for  selecting  the  unique  best.  The  procedures  include  preliminary  tests  which 
allow  the  experimenter  to  have  an  option  to  not  select  if  the  statistical  evidence  is  not 
significant.  Two  probabilities,  the  probability  to  make  a  selection  and  the  probabil¬ 
ity  of  a  correct  selection,  are  controlled  by  these  selection  procedures.  Comparisons 
between  the  proposed  selection  procedures  and  certain  earlier  existing  procedures 
are  also  made.  The  results  show  the  superiority  of  the  proposed  selection  procedures 
in  terms  of  the  required  sample  size.  ^  ^ 


1.  INTRODUCTION 


The  problem  of  selecting  the  best  population  from  among  k{>  2)  populations  has 
been  studied  extensively.  A  lot  of  selection  procedures  have  been  derived  for  different 
selection  goals  by  several  authors.  Among  them,  Bechhofer  (1954)  introduced  the 
indifference  zone  approach  for  selecting  the  normal  population  with  the  largest  mean. 
In  his  approach,  the  determination  of  the  sample  size  is  heavily  dependent  of  the 
indifference  zone  assumption.  Also,  the  probability  of  a  correct  selection  depends 
on  the  unknown  parameters  and  is  analogous  to  the  power  of  a  test.  However,  a 
probability  that  is  analogous  to  the  probability  of  type-I  error  of  a  test  was  not 
taken  into  consideration.  Bechhofer’s  procedure  forces  the  experimenter  to  make 
a  selection,  and  often  that  procedure  is  not  used  in  practiced  application  because 
of  the  lack  of  a  statistical  test  for  the  homogeneity  of  the  parameters  (see  Simon 
(1977)). 

Based  on  the  preceding  reasoning.  Bishop  and  Pirie  (1979)  introduced  a  selection 
procedure  in  which  a  test  of  homogeneity  was  conducted.  The  procedure  allows  the 
experimenter  to  have  the  option  not  to  make  a  selection  if  the  statistical  evidence  is 
not  significant.  Later,  Chen  (1985)  proposed  modified  selection  procedures  for  the 
problem  of  selecting  the  best  normal  population.  He  considered  a  preliminary  test 
based  on  the  sampled  spacing  between  the  largest  and  the  second  largest  ordered 
statistics.  If  the  statistical  evidence  of  the  preliminary  test  is  not  significant,  the 
experimenter  decides  not  to  make  a  selection.  Otherwise,  he  selects  the  population 
yielding  the  largest  sample  mean  value  as  the  best  population.  The  sample  size 
is  determined  to  control  both  the  probability  of  type-I  error  for  the  preliminary 
test  and  the  probability  of  a  correct  selection.  Analogous  to  Chen  (1985),  Chen 
and  Mithongtae  (1986)  proposed  selection  procedures  for  two- parameter  exponential 
distribution  models.  However,  their  procedures  are  conservative  in  the  sense  that 
the  determined  sample  size  is  always  larger  than  the  required  minimum  sample 
size.  Also,  their  procedures  can  not  be  applied  for  a  case  where  the  common  scale 
parameter  is  unknown. 

In  this  paper,  along  the  line  of  Chen  (1985)  and  Chen  and  Mithongtae  (1986), 


we  study  the  problem  of  selecting  the  best  population  from  among  k{>  2)  two- 
parameter  exponential  populations.  Selection  procedures  are  derived  according  to 
whether  the  common  scale  parameter  is  known  or  unknown.  Exact  sample  sizes  are 
determined  to  control  both  the  probability  of  type-I  error  and  the  probability  of  a 
correct  selection  for  the  cases  considered  by  Chen  and  Mithongtae  (1986).  When 
the  scale  parameter  is  unknown,  a  two-stage  selection  procedure  is  also  proposed. 
The  two-stage  selection  procedure  covers  the  case  where  the  procedures  of  Chen 
and  Mithongtae  (1986)  can  not  be  used.  It  should  be  pointed  out  that  the  proposed 
procedures  can  be  implemented  by  using  certain  existing  tables. 

2.  FORMULATION  OF  THE  PROBLEM 

Let  denote  k{>  2)  independent  two-parameter  exponential  distribu¬ 
tions  with  guaranteed  life  times  respectively,  and  a  common  (known  or 

unknown)  standard  deviation  0.  Let  //[ij  <  . . .  <  /i[*j  denote  the  ordered  values  of 
/xi , . . . ,  /Zfc  in  the  parameter  space  ft  where 

ft  =  {(^,0)l/£  =  (/ii,...,Mfc),-oo  <  fii<  oo,  e>  0}. 

We  partition  the  parameter  space  into  the  following  three  subspaces: 

the  preference  zone;  Sl{PZ)  =  {(£,^)  €  ft|  >  6,  S  >  0}, 

the  nonselection  zone:  il{N Z)  =  {{l£,0)  G  ft|^[jt_i]  = 

and 

the  indifference  zone:  il{IZ)  =  ft  -  ft(FZ)  -  il{NZ), 
where  ^  is  a  known  positive  constant. 

Denote  the  event  of  a  correct  selection  by  CS  and  the  event  of  selection  by  S. 
The  goal  is  to  develop  selection  procedures  R  to  select  a  single  best  population  with 
a  minimum  sample  size  from  each  population  such  that  the  following  probability 
requirements  are  satisfied: 

P^^,e){S\R)  <  a  for  all  (^,0)  G  ft(AZ)  (2.1) 

and 

Pi^,0)iCS\R)>  P*  for  ail  ill,  0)  G  ft(PZ)  (2.2) 


The  selection  procedure  R  depends  on  whether  the  common  standard  deviation 
is  known  or  unknown. 

3.  SELECTION  PROCEDURE  FOR  STANDARD 
DEVIATION  KNOWN  CASE 

Let  Xii,  • .  • ,  be  a  sample  of  size  n  arising  from  population  TTj  and  let  Yi  = 
^min  Xij,  i  —  Also,  let  Vji]  <  ...  <  yjk)  denote  the  ordered  statistics  of 

Pi, . . . , yjfc.  When  $  is  known,  we  propose  a  selection  procedure  as  follows: 

Ri:  select  the  population  yielding  Vjfc]  as  the  best  population  if  yj/t)  -  ^[^-1]  > 
Xd/n;  otherwise,  do  not  make  a  selection,  where  A  and  n  are  chosen  to  satisfy  the 
probability  requirements  (2.1)  and  (2.2). 

For  the  rule  Ri,  we  need  to  investigate  the  supremum  of  )  for  (£,  9)  6 

9.{NZ)  and  the  infimum  of  P(^,j)(CS|Ri)  for  (£,^)  6  fl{PZ). 

Firstly,  let  Afj,  t  =  1, . . . ,  fc  be  independently  distributed  as  two-parameter  expo¬ 
nential  distribution  with  guaranteed  life  time  1, . . . ,  A:  and  common  standard 

deviation  1  respectively.  Let  ATjii  <  ...  <  Jfi/fc]  denote  the  ordered  statistics  of 
Xi, . ,  .,Xk,  A*!!)  <  ^  M(fcl  be  the  ordered  guaranteed  life  times,  and  X(,)  be  the 

random  variable  associated  with  parameter  /ijjj,  t  =  1, . . . ,  fc,  respectively.  The.\  we 
have  the  following  Lemma: 

Lemma  3.1.  For  any  positive  constant  c,  >  c}  is  nonincreasing  in 

Uli],  i=  l,...,fc-2. 


Proof.  Let  F  be  the  cdf  and  /  be  the  pdf  of  the  standard  exponential  random 
variable.  By  symmetry  we  may  assume  that  fii  <  fi2  <  . . .  <  fik-  For  any  c  >  0,  we 
have 


where 


^  _  J  t^k  1, 

'  "  (  max(0,//*:-i +c),  i  =  k. 


For£=l,2,...,k-2, 

^PufXm  -  x,,.,,  >  cj 


c)f(x  +  in-  Hi-  c)f{x)dx 


0OO  ** 

=  Y^  Yl  F{x + fit  -  - 

,=1  j=l 

i^t 

k^i  »oo 

-  ^  /  JJ  /’(a:  +  -  Hj  -  c}fix  +  /r,  -  fig-  c)f{x)dx 

t=i  j=i 

J^i.t 

1^00 

-  /  Yl  Hx  +  ftk  -  -  c)fix  +  fik  -He-  c)f(x)dx 

Jdk  j_i 


i^t 

=  h  +  h  -  hi 


(3.1) 


where 


h 


k-l  ,c 

-u. 

i*e 


( 


ji  F(a;  + /if  - /ij  +  c) 


\m.t 


[f{x  +  H(-  Hi  -c)/(x) 


y 


and 


-f 


-r 


k-1 


-  f{x  -  c)/(x  +  /if  -  fdi)]dx, 
\ 


n  +  fte-ftj-c) 


k-l 


f{x  +  /if  -  /ifc  -  c)/(i)di, 


]][  F{x  +  /ifc  -  Hj  -  c) 


\i/< 


f{x  +  Hk-  fit-  c)f{x)dx. 


Ii  =  0,  since  /(x  +  /if  -/ii  -c)/(x)  =  /(x-c)/(x  +  /if-/ii)-  For  I2-I3,  we  consider 
the  following  two  cases: 

Case  L  Hk-i  -  Hk  +  c>0. 


After  changing  variables, 


I2  -  h 


-f 

j  Uk- 


fc-1 


J|  F{x  -\rtH-  fij  -  c)  [f(x  c)f{x)  -  /{x  -  c) 


fUk-Mt+c 
JMk-l-tkt  +  C 


f{x  +  nt-  fj’k)]dx 
/.  .  \ 


fc-i 


n  + M/  -  Mi  -  c) 


f{x  -  c)f{x  He  -  fik)dx.  (3.2) 


) 


The  expression  in  the  bracket  of  the  first  term  in  (3.2)  is  zero  and  the  second  term 
is  nonnegative,  hence  the  derivative  (3.1)  is  nonpositive. 


Case  2.  Hk-\  -  Mfc  +  c  <  0. 


After  changing  variables,  analogous  to  the  argument  of  Case  1,  we  have 

h-h 

n  -  Mi  - 


-/ 

J Hk  -M<+C 


y 


[f(x-^Ht-  ^^k-c)f(x) 

-  /(*  -  c)/(x  +  Me  -  f^k)]dx 

f{x  -  c)f{x  +  He  -  fJ'k)dx 


( 

\ 

fttk -/it  +  C 

k-\ 

/ 

JJ  F(x  +  He  - 

-  Mi  -  c) 

-'M* -Ui 

/ 

<  0. 

This  completes  the  proof  of  Lemma. 


□ 


From  Lemma  3.1.,  we  have  the  following  corollary: 


Corollary  3.2.  PjA^\k]  -  >  c}  ^  P{\^i  -  ^2!  >  c}  for  all  ^  G  il{N Z),  where 

Zi  and  Z2  are  independent  having  standard  exponential  distribution.  Furthermore, 

sup  PAX\k]  -  -X'(fc_i)  >  c}  =  2H{-c), 
n(NZ) 

where  H(t)  denotes  the  cdf  of  Z\  —  Z^. 


Note  that 


1  - 


2®  ) 


for  <  >  0 
for  t  <  0. 


(3.3) 


Since  n{Yi  -  fti)/6  is  distributed  as  a  standard  exponential,  we  have  the  following 
result: 

Iheorem  3.3.  sup 

OiNZ)  - 

P-£Q9f.  For  (£,0)  G  n(iVZ), 

=  P(ii,e){y[k]  -  Vifc-i]  >  XO/n) 

-  - ^  >  A} 

^  ~  ■^2!  >  A},  by  Corollary  3.2. 

=  2Hi-X)  =  e-\  □ 


R^niftrk  3-li  In  order  to  satisfy  the  probability  requirement  (2.1),  we  may  let 
e~^  =  Q.  Thus 

A  =  -In  a.  (3.4) 

We  next  evaluate  the  infimum  of  P^^^ey{CS\Rl)  over  il(PZ). 

Thwrglll  3i4i  The  infimum  of  over  U{PZ)  occurs  at  the  configura¬ 
tion  ^[1]  =  . . .  =  =  /iifc]  -  60  and 

n{pz)  =  /"(I  -  (3.5) 

where  t  =  max(0,  A  -  r)  and  r  =  n6.  Further,  if  r  >  A,  then 

-(1  -c-^"-"))*)//!:.  (3.6) 


Proof.  See  Chen  and  Mithongtae  (1986). 


□ 


Remark  3.2. 

(a)  If  r  <  A,  then  the  right  hand  side  of  (3.5)  is  which  is  less  than  or  equal 

to  l/k. 

(b)  Since  r  =  nS,  for  n  large  enough,  we  have  r  >  A.  In  order  to  satisfy  the 
probability  requirement  (2.2),  from  (3.6),  we  may  let 

e^-^(l-(l-e-<^-^))'')/it  =  P*.  (3.7) 


(c)  Let  V  —  ,  then  (3.7)  can  be  rewritten  as 

t,-i(l_(l_„)*)  =  ;fcP*.  (3.8) 


Tables  for  the  v  values  when  P"  =  0.90,  0.95,  0.975  and  0.99;  k  =  2(1)  15  are 
available  from  Raghavachari  and  Starr  (1970). 

(d)  For  a  £  (0, 1),  P*  £  (1/fc,  1)  and  k,  let  A  =  -£n  a.  Also  let  v  be  the  solution  of 
the  equation  (3.8).  Then  r  =  —in{av).  Let  <  x  >  denote  the  smallest  integer 
not  less  than  x,  then 

n  =<  —£n{av)l6  >  (3.9) 

is  the  smallest  sample  size  satisfying  the  probability  requirements  (2.1)  and  (2.2). 

(e)  In  order  to  compare  our  procedure  with  Chen  and  Mithongtae’s  procedure,  we  let 
n*,  A*,r*  and  6"  denote  the  corresponding  notations  of  Chen  and  Mithongtae’s 
procedure. 

Then 

A*  =  A  +  u,  r*  =  n'6,  6*  =  60,  t  =  n6 


and 


k-1 


where  «  =  i  ^  1/(A:  —  i -I-  1). 


»=i 


Hence 


n*  =  n  +  {in  u)f6. 


Thus,  n*  >  n  as  fc  >  3. 

(f)  The  Table  1  of  Chen  and  Mithongtae  (1986)  can  be  used  to  implement  our 
procedure.  The  minimum  sample  size  required  is  n  =<  -£n(av)j6  >,  where  v  = 
and  A*  and  r*  can  be  found  from  Table  1  of  Chen  and  Mithongtae 

(1986). 

Example.  For  k  =  5,  ff  =  2,  6“  =  0.5,  a  =  0.05  and  P*  =  0.95,  then  n’  =  29  and 
n  =  27. 

4.  SELECTION  PROCEDURE  FOR  STANDARD 
DEVIATION  UNKNOWN  CASE 

Under  the  formulation  of  Q.{PZ)  in  Section  2,  we  have  a  single  stage  selection 
procedure  when  the  standard  deviation  is  unknown.  We  take  a  sample  of  size  n,  say 
ATjj,  j  =  l,...,n  from  each  population  tt,,  i  =  l,...,fc.  Let  y,  be  the  first  order 
statistic.  The  minimum  variance  unbiased  estimator  of  B  is 

k  n 

<=i  j=i 


where  v  =  k{n  -  1).  Let  Fii]  <  ...  <  K(fc]  denote  the  ordered  values  of  y,’s.  We 
propose  a  selection  procedure  as  follows; 

Ri:  select  the  population  yielding  Fj*]  as  the  best  population  if  >^[*1  -  > 

'yO/n;  otherwise,  do  not  make  a  selection,  where  7  and  n  are  chosen  to  satisfy  the 
probability  requirements  (2.1)  and  (2.2). 

Note  that  Zi  =  n(Y^i)  -  fi\{\)/B,  i  =  l,...,k  and  B  are  mutually  independent. 
Zi  ~  standard  exponential  and  Q  =  B/B  has  a  gamma  distribution  with  parameters 
V  and  V.  We  use  G{x)  to  denote  the  cdf  of  Q.  Then  g{x)  =  G'{x)  = 

By  applying  Corollary  3.2,  we  have  the  following  result: 


Theorem  4.1. 


(4.1) 


sup 

n(/vz) 


Proof. 


Now 


Hence 


-  P{yik]-yik-n  >  l9ln) 

=  E{P(^^^e){yik]  -  yifc-il  >  iQ&I'nlQ})- 

P{ji,e){^lk]  ~  ^lk-i|  >  7Q^/«IQ  =  a;} 

=  >  IxO/n] 

<  2H{~‘yx),  by  Corollary  3.2. 


sup  P(f,,e)iS\R2} 
n(NZ)  “ 


-i: 

-L 


e-'*^dG{x) 


OO  ^v^v-1 


=  ( 


7  +  U 


r(u) 

y. 


□ 


By  setting  sup  9)(51iZ2)  equal  to  q,  we  obtain 
n(NZ) 

7  =  -  1).  (4.2) 

Now  we  evaluate  the  infimum  of  P(^,j)(C5|/?2)  for  (/£,  ^)  €  fl(PZ). 

Theorem  4.2.  The  infimum  of  P(^,j)(C5|i22)  over  il{PZ)  occurs  at  the  configura¬ 
tion  ^(1)  =  . . .  =  =  (i[k\  -  dO  and 

/nS/7  1 

JillAliMCS\R2)=  / 

*•(  " X)  Jq  K 

rOO  1 

+  /  14.31 


Remark  4.1.  Although  Theorem  4.2.  is  similar  to  Theorem  4  of  Chen  and  M  ithong- 
tae  (1986),  the  concept  is  different.  Their  procedure  should  specify  the  value  of  b*  19 
in  advance.  However,  it  seems  hard  or  not  possible  since  the  value  of  the  parameter 
6  is  unknown. 

In  order  to  satisfy  the  probability  lequirement  (2.2),  we  set 

rns/i  ,  roo  , 

/  -i€"«-^=^(l-(l-c^^-"^)*)dG(i)+  /  je^^-^^dG{x)  =  P-.  (4.4) 

Jo  *  JnSf-y  * 

The  numerical  methods  of  Chen  and  Mithongtae  (1986)  can  be  applied  for  solving 

equation  (4.4). 

5.  A  TWO-STAGE  SELECTION  PROCEDURE  WITH 
UNKNOWN  STANDARD  DEVIATION 

In  the  following,  we  consider  the  case  where  the  parameter  space  is  partitioned 
by  the  same  way  as  that  of  Chen  and  Mithongtae  (1986).  That  is, 

(l{PZ)  =  {(£,^)|M[fci  -  >  ^,  ^  >  0), 


and 

iliNZ)  =  {(£,0)l/i(*_i)  =  ^[*]}. 

When  0  is  known,  the  situation  is  equivalent  to  the  one  discussed  in  Section  3. 
However,  when  9  is  unknown,  a  single  stage  selection  procedure  proposed  by  Chen 
and  Mithongtae  (1986)  can  not  be  implemented  since  it  seems  not  possible  to  specify 
the  value  of  the  ratio  6/9  in  advance.  Analogous  to  that  of  Desu,  Narula  and 
Villarreal  (1977),  a  two-stage  selection  procedure  is  proposed  as  follows: 

(i)  Take  an  initial  sample  vif  size  no  from  each  of  the  k  populations,  say  A,i,..., 

.  fc  no 

Xing,  i  =  Let  Yi{no)  =  min  Aij,  t  =  l,...,fc,  and  0o  =  E 

Yi{no))/vo,  where  vo  =  k{no  -  1). 

(ii)  compute  rj  =  Mq/S,  where  h  >  70  is  determined  so  that 

r  /°°(i  _  =  P\ 

Jo  Jo 


(5.1) 


7o  =  -  1}  and  g(x)  is  the  pdf  of  Q  =  0o/^-  Let  N  =  m<Lx{no,  <??>}, 

where  <  x  >  denotes  the  smallest  integer  not  less  than  x. 


(iii)  If  necessary,  take  N  —  no  additional  observations  from  each  of  the  k  populations 
and  compute 

YJN)  =  min  X,,,  i  =  l,...,k. 

(iv)  Select  the  population  yielding  yj;t|(7V)  as  the  best  population  if  yiA;](X)  —  VjAr-ij 
{N)>  'yoOo/N;  otherwise,  do  not  make  a  selection. 

Let  Zi  =  N{Yi(N)  —  i  =  It  is  well-known  that,  given  N  = 

n,  Zi,  i  =  1, . . .  ,k,  are  i.i.d.  having  standard  exponential  distribution. 

Now,  we  evaluate  the  supremum  of  P(^^e)(Sli?3)  for  (^,0)  G  Cl(NZ). 

Theorem  5.1.  For  any  (^,0)  G  Q(NZ)  and  any  positive  constant  70,  P(^,e){y|fc](7V) 
-y(k-i](X)  >  to^o/A"}  <  where  vo  =  k(no  -  1). 

Eraaf,  Let  An  =  {N  =  n},  then  =  {Q  <  ^}  and  An  =  <Q  < 

n  >  no. 


^(ii.«){51fcl(iV)  -  llfc-il(A^)  >  7o^o/iV} 


>  7oQ 


nssno 


00  - 

<  ^  /  e~'^°^g{x)dx,  by  corollary  3.2. 

n=no 

=  /  e~''^’^g{x)dx 
Jo 

V70  +  V0/ 


Remark  5.1.  Note  that  sup  which  is  independent  of  iV. 

ft(ffZ) 

In  order  to  satisfy  the  probability  requirement  (2.1),  we  may  set  )"°  =  O’  and 


hence 


7o  =  no{a'’/*'o  -  1}.  (5.2) 

Next,  we  evaluate  the  infimum  of  the  probability  of  a  correct  selection  for  (£,  0)  € 
fi(PZ). 

Theorem  5.2.  The  infimum  of  P(^^s)(C5|P3)  over  Cl{PZ)  occurs  at  the  configura¬ 
tion  /f[ij  =  . . .  =  =  ^l[k]  -  6  and 

inf  P(^,e)(C51P3)  >  /  /  (5.3) 

{l(.PZ)  -  Jo  Jo 

PlQOf,  P(^,tf)(C5|J?3) 

=  /'(M){y(*)(^)  -  Yik^,]{N)  >  70  w  Y^k){N)  =  Y^k]{N)} 

=  P{ti,e){Y{j){N)  <  Y^k){N)  -  -toffolN,  j  =  1,. .  .,fc  -  1} 

_  n  (^(yum-rn)  .  N{Y^k){N)-^Hk])  ,  N{fHk]-m) 

0  ^  9  0 

>  P^Zj<Zk  +  ^  -  7o<?>  j  =  1,...,*:-  l| 

N6 

>  P{Zj  <  Zk  +  hQ  -  7oQ,  j  =  -  1},  since  -j-  >  hQ 


If  ft  <  ^0)  then  the  right  hand  side  of  (5.3)  is 


noo 

(1  -  ^ e~‘'d(g{x)dx 

,70-/i)x 

noo 

=  -f 

kJo 


p(A-7o)ii!o 


vX° 


e-'’°^dx 


r(wo) 

=  .y 

k  \7o  -  h  +  voj 
<  since  7o  -  h  +  vo  >  vq. 

Hence  we  should  choose  ft  such  that  ft  >  70. 

In  order  to  satisfy  the  probability  requirement  (2.2),  we  may  set 

^00  ^00 

I  I  (1  _  ^  e~^d(^g{x)dx  =  P’ . 

Jo  Jo 


(5.4) 


Remark  5.3.  Table  of  Desu,  Narula  and  Villarreal  (1977)  can  be  used  to  implement 
the  proposed  two-stage  procedure  R3.  Let  d  be  the  value  obtained  from  their  table, 
then  ft  =  7o  +  2vod. 
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